gr-qc/9911062 



UCVFC-DF/16-99 



Analytical Description of Voids in 
Majumdar-Papapetrou Spacetimes 

Victor VarelaQ 

Department of Mathematics and Statistics, 
The University of Edinburgh, 
Edinburgh, Scotland 

Escuela de Fisica, Facultad de Ciencias, 
Universidad Central de Venezuela, 
Caracas, Venezuela Q 
E-mail:varela(§belka . ciens . ucv . ve 

Departamento de Fisica, 
Universidad Simon Bolivar, 
Caracas, Venezuela 

Abstract 

We discuss new Majumdar-Papapetrou solutions for the 3+1 Einstein-Maxwell 
equations, with charged dust acting as the external source of the fields. The solutions 
satisfy non-linear potential equations which are related to well-known wave equations of 
soliton physics. Although the matter distributions are not localised, they present 
central structures which may be identified with voids. 
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1 Introduction 



We consider solutions for the Einstein-Maxwell (EM) equations with charged dust 
acting as the external source of the fields. Our basic equations read 

G'i = SttTI:, (1) 

Ff"".,, = AttJ'', (2) 

where and F^^'^ denote the Einstein and Maxwell tensors, and the total energy- 
momentum tensor is given by 

= Ej^ + pu^'u,. (3) 

Here E^ is the Maxwell energy-momentum tensor, and the matter term corresponds 
to dust with energy density p and four- velocity u^. The four-current is defined by the 
expression 

= au^, (4) 

were a is the charge density. 

We assume that the fluid is static and use the conforstatic metric 

ds'^ = -y'^ dt^ j^^hijdx'dx\ (5) 

where the background metric hij and V depend only on the space-like coordinates x^, 
x^, x^. The electrostatic forms of A^^ and are given by 

= Aoix')5'^, (6) 

with i = 1, 2, 3. 

Under these conditions, Eq. (0) contains only one non-trivial equation: 



(8) 



where h and /i*-' are the determinant and the inverse of hij, respectively. 
The trace of the Einstein equations is 

R = -SttT, (9) 

where R denotes the Ricci scalar and T = Tj^. We use the decomposition 



R = V 



Rh + 2Vi In y - 2di In Yd' In V 



(10) 



Here Rh is the Ricci scalar associated to hij, and V\ is the three-dimensional Laplacian 
operator constructed with the same metric. We assume a flat background space, with 
Rh = 0. Therefore, combining Eqs. (|9|) and ( p!o| ) we obtain 



^^^ =17^- (11) 
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Following the Majumdar-Papapetrou (MP) procedure, [|l|, ^ we assume that 

^0 = aV, (12) 
where a = ±1. As a consequence, the Maxwell equation (^) takes the form 



which is clearly the same as Eq. ( [TlD whenever the condition 

T = -aJ^V (14) 
holds. This equation can be combined with j'' = y to obtain the alternative expression 

a = -aT. (15) 
bmce 1 = —p for dust, Eqs. (0) and (ID can be finally expressed as 

V^A + iirpX^ = 0, (16) 

a = ap, (17) 

where A = p. Due to Eqs. (^) and ([121), only one Einstein equation is not trivially 
satisfied. ^ Therefore, solving Eq. ( [iq) is sufficient for finding a solution of the EM 
equations. 

If we identify our fiat background space with the Euclidean, three dimensional space 
and assume p = 0, then Eq. (|T^) reduces to the usual Laplace equation V^A = and 
the electrovac, multi-black hole solution follows straightforwardly. Assuming spherical 
symmetry, and using spherical coordinates, we find 

A = l + -. (18) 
r 

In the far-asymptotic region, the behaviour of this solution is approximately given by 

y«l--, 500^-1 + ^, Ao«±(l--). (19) 
r ' r 

The corresponding expression for the electric field is 

E^S^, (20) 

where 



q = ±m. (21) 

Equation (^) implies that the invariant area of any 2-sphere surrounding the origin 
is given by y^p^- Therefore, the set r = 0, t = constant has a non-zero invariant area 

given by ^Trm?. In fact, a simple coordinate transform shows that the null hypersurface 
r = is the horizon of the extremal Reissner-Nordstrom solution. Also, if we define 
the new radial coordinate f = —r and perform the standard analysis, ||5|] then we find 
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that this horizon encloses a point-like, essential singularity placed at f = m. In fact, 

the invariant area vanishes and the scalar J = F^yF^^'^ = (^) blows up at that 
point. 

Equations (|l^) and ( p^ were originally discussed by Das in his study of equi- 
librium configurations of self-gravitating, charged dust. More recently, Giirses has 
considered non-electrovac solutions when Eq. (^) is linear. This situation corresponds 
to his choice p = for constant h. In this case, Eq. (|l^) admits the particular so- 
lution A = ^^^^ where a is an integration constant. The oscillatory behaviour of 
this solution implies a geometry with a complicated radial dependence. In fact, the 
invariant area vanishes for a discrete, infinite set of values of r, and the Ricci scalar 
R = , blows up wherever the invariant area vanishes, except for r = 0. Other 

sin or ^ 

solutions with oscillatory behaviour have been considered by Balakrishna and Wali,[0] 
Braden and Varela,|8| and Ida.|Q] In Section 2 we exploit the general non-linearity of 
Eq. (|l|) to obtain new solutions which are free of oscillatory singularities and allow 
asymptotically flat behaviour. 



2 The non-linear models 

The non-linear potential equation (p!6[) takes the spherically symmetric form 



(fx 2dX n 

+ - +4^pA3 = 0. 22 
ar^ r ar 



Using the new radial coordinate r = the same differential equation can be written 
as 



If p and A satisfy the condition 



5^ sin A 

then ( p^ ) finally reduces to the -sine-Gordon equation |^ 

^+62sinA = 0, (25) 

which has the solutions 

A^ (r) = 2 arcsin [tanh (±6t + c)] + 2mT, (26) 

where n is an arbitrary integer, c is an integration constant, and b is assumed to be 
positive. We consider only the case n = 0. In terms of the original radial coordinate, 
these solutions read 

(0 = r ■ / . XT • (27) 



2 arcsin 



tanh [±- + c 
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We observe that V^{0)'^ is finite, so the invariant area vanishes for r = 0. Therefore, 
the set r = 0,t = constant is point-hke with respect to both solutions. Let us deal with 
first. A preliminary numerical study of the invariants J, R, R^'^Rap, R°'^'^^ Ra/s^yS 
suggests that these quantities are bounded for non-negative r, whenever c is positive. 
If we choose 
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1 + sin(l/2) 



.l-sin(l/2)J' 

then the far-asymptotic behaviour of this solution is given by Eqs. (|T9|), (pO|), (21) with 
m = 26cos(l/2). Therefore, is asymptotically flat, exactly as the MP electrovac 
solution. 

The positive definite energy density given by Eq. ( p^ ) corresponds to a non-localised 
matter (and charge) distribution. However, p is negligible for x = | ^ 0.2. For very 
small X the dimensionless expressions of p and R"'^^^ Rap-yS are approximately given by 

-c -i 

P(^) ^ ^ ^> (29) 
vr* x^ 

6 — - 

i?"^^'i?a;375W«3(^^) e-^'^. (30) 

These results imply a very fast decrease of p{x) and R"^'^^ Rai3js{x) when x — > 0"*", and 
suggest the existence of a void in the innermost region of the asymptotically flat object 
constructed with V~^. Nevertheless, this interpretation cannot be complete without a 
better understanding of the point-like set r = 0, t = constant. A closer look at the 
singularity contained in this solution is also necessary. We interpret r = 0, t = constant 
as the center of symmetry and observe that the above mentioned invariants are bounded 
at this point. However, the coordinate transform f = — r reveals the existence of a 
point-like, essential singularity at f = ^. In fact, the invariant J blows up at this 
point. The use of alone may imply the division of the manifold into connected 
parts, separated by the point-like singularity placed at r = — |. However, a very 
different situation comes out when we restrict to positive values of r and describe 
the geometry for r < with the second solution V~ . Then, a smooth (at least C^) 
matching of and occurs at r = and the arising asymptotically flat spacetime 
seems to be connected and singularity free, with an almost empty region near the center 
of symmetry. Thus, the joint use of and V~ provides a simpler description of a MP 
void. The study of the global structure of these solutions is left as an open problem, 
which provides motivation for further research work. 

Finally, we point out that other exact, non-linear solutions for this theory can be 
found if we impose different relationships between p and A. For example, the choice 

ft^r^sinA 

leads to the sine-Gordon equation 

^ = 52 sin A. (32) 
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It has the well-known solutions 

A±(t) =4arctane(±''"+'^). (33) 

If we choose d = In [tan(l/4)], then both solutions have asymptotically flat behaviour. 
Another example is 

p=^(A-A3). (34) 
In this case the geometry is determined by the A(/>^ equation 



+ 62 (a - A^) = (35) 



A±(t) = tanh ( ±4=^ + / ) • (36) 



which admits the solutions 

\'^(t'\ — fnnV, ( 

The relationship between the 3+1 EM theory and the equations of 1+1 soliton 
physics deserves a more detailed examination. Possible extensions of this work involve 
the analysis of dust models for which A(t) is a solution of the KdV equation, and the 
study of the non-linear potential equations arising in higher dimensions.! 
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